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Abstract 



We introduce a toy model for interacting populations connected by mutations 
and limited by a shared resource. We study the presence of Eigen's error 
threshold and mutational meltdown. The phase diagram of the system shows 
that the extinction of the whole population due to mutational meltdown can 
occur well before an eventual error threshold transition. 
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The evolutionary process is due to the balance between two opposite forces: random 
mutation and natural selection. Mutations tend to increase the genetic diversity, while 
selection discriminates the individuals more apt to survive. A simple model for evolution of 
self- reproducing molecules was proposed by Eigen |j] . The inheritable characters (genotype) 
are modeled by a binary string g = (gi, g 2 , . . . ,gi), 9% £ {0, 1}, of fixed length L (haploid 
individuals). The genotypic space is thus a binary hypercube of L dimensions. One can 
think of the zeros as good genes (or basis) and the ones as bad genes. We shall denote the 
set of individuals sharing the same genotype as a strain. 

The selection is modeled by the concept of fitness, which can be defined as the survival 
probability or birth rate of a strain in the limit of zero mutations and vanishing population 
(in order to avoid overcrowding effects). The selection acts on the external characters of an 
individual, i.e. its phenotype m, which can be considered a (generally degenerated) function 
of the genotype m(g). The simplest phenotype (used also in Eigen's model) is a decreasing 
function of the fraction of bad genes in the genotype: m(g) = (1/L) J2i 9% ■ The mutations, 
occurring with probability /i, reverse the value of a single bit (point mutations) and thus 
couple different strains. 

The classification of strains into species is based both on the phenotypic traits and on 
the genotypic information: for strains to belong to the same species, they must be connected 
by mutations and form an isolated cluster in the phenotypic space. This last requirement is 
guaranteed, even in an almost flat fitness landscape, by competition ||, which "eliminates" 
the less adapted strains phenotypically near to a fittest one. 

For a sharp peak landscape the genotype g = (0,0,...) has higher fitness a, and all 
other genotypes have the same lower fitness (3: the fitness landscape is almost flat. Two 
different regimes are possible: in one phase the asymptotic limit is a bell shaped distribution 
in phenotypic space (quasispecies), centered around the master sequence; in the other the 
most common phenotype is no more the fittest one. The shape of the distribution in this 
second phase is dominated by combinatorial factors, approaching a Gaussian in the limit of 
very large mutation rates, with almost no difference among strains. 

The transition between these two regimes, triggered by the mutation rate /i or by the 
length of the genome L, is called error threshold. While the original results have been 
obtained for an infinite population, this transition, disregarding the effects of fluctuations, 
depends slightly on the size of population ||^[]. The error transition has the character of a 
true phase transition (first order) for an infinite genome length L |J. 

A related effect, effective in small populations, is called Muller's ratchet or stochastic 
escape [0,§: in absence of back mutations there is a finite probability that the master 
sequence will be lost due to fluctuations. Since it relies on a random process, the average 
escape time is however of order of the exponential of the size of the strain , and thus this 
effect is relevant only for very small populations. In practice, this is the ultimate extinction 
mechanism, effective when another cause has reduced enough the population. The error 
threshold, on the other hand, does not depend on the size of the total population. In the 
following we shall neglect to consider the influence of fluctuations. 
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When the total population size is not kept fixed, the presence of unfavorable mutations 
reduces the average fitness, and this can lead to the extinction of the whole population, a 
phenomenon which is called mutational meltdown JTOHl^]. The presence of this effect in 



the Eigen model has been shown numerically in Ref. JT3[ . 

In this work we propose a simple minimal model which is able to exhibit both the 
error threshold and the mutational meltdown. Due to the form of the fitness function, 
the dynamics of the population is fundamentally determined by the fittest strains. Let us 
indicate with X the number of individuals sharing the master sequence, with Y the number 
of individuals whose genotype has m = 1, and with Z all others individuals. We assume 
also non-overlapping generations, so that we can consider a discrete time dynamics with 
unit equal to the generation time. 

During reproduction, individuals from strain X can mutate, contributing to Y, and Y 
to Z, with mutation rate /iL. We disregard the possibility of back mutations from Z to Y 
and from Y to X. This last assumption is equivalent to the limit of a large genome, which 
is the case for existing organisms. We also introduce the reproduction accuracy q, defined in 
term of fiL as q = 1 — fiL. Due to the assumption of large L, the multiplicity of mutations 
from m — 1 to m > 1 (L — 1) is almost the same of that from m = to m — 1 (L). 

We shall assume a finite (and constant) carrying capacity K of the environment, assuming 
that the effective reproduction rate of a population is proportional to 1 — N/K, where 
N = X + Y + Z is the total number of individuals. The evolution equation of the population 
is 
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(1) 



where the prime denote quantities at following time step. Obviously all population dies if 
a < 1. 

By introducing the normalized distribution x = X/N, y = Y/N and z = Z/N (x+y+z = 
1) we have 
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and 
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X' = X { 1-— J (ax + P(l-x)). 
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The steady state of Eq. (0) is given by x' 
fixed points 



and 
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The fixed point A is always unstable and B always stable for a > 1. In this second case the 
asymptotic population size is 
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We characterize the error threshold by the fact that the fittest phenotype is no longer 
the most common one, its phase boundary is given by x = y. Evaluating this condition at 
the fixed point B we get from Eqs. (|5|) 
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The mutational meltdown corresponds to the vanishing of the total population N 
while not changing its distribution. Thus, from Eq. (0) we have 
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(7) 
= 0, 

(8) 



Our definition of the error threshold transition needs some remarks: in the original 
work the error threshold is located at the maximum mean Hamming distance, which 
corresponds to the maximum spread of population. In the limit of very large genomes these 
two definitions agree, since the transition becomes very sharp ||. For finite populations, a 
possible transition can be located in correspondence of the vanishing of the probability of 
finding the master sequence in a way similar to Muller's ratchet. 

However, since the master sequence has always an higher effective fitness than other 
strains, this implies that for (3 < 1 the error threshold is only a transient phenomenon before 
extinction, unless the population is artificially kept finite, as in Ref. Pj. This mutation- 



induced extinction has been investigated numerically by Malarz and Tiggemann in Ref. |L3 

The condition for the disappearance of the master sequence is x'/x < 1 in the limit of 
vanishing x. Denoting arbitrarily this transition as Muller's ratchet's one, we obtain from 
the first of Eqs. (H) 
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This transition coincides with the mutational meltdown for [3=1. 

In Figure 1 we plot the phase diagram of the three transitions for some different values 
of (3. We note that the error threshold depends only on the ratio a/ (3 (i.e. on the relative 
fitness of strains); on the other hand, the mutational meltdown threshold involves the whole 
population, so it depends on the absolute value of the fitness a. For (3 < 1 the mutational 
meltdown transition crosses the error threshold, so one can observe both extinction of the 
quasispecies distribution or the error threshold in a stable population, according with the 
parameters. 

The physical relevant conditions for the extinction of the master sequence is (3 < 1, for 
which also the Y strains vanish; otherwise back mutations cannot be neglected. Our toy 
model can approximate the M2 model by Malarz and Tiggemann, in which only one bit 
is mutated with probability P mu t, and thus q = 1 — P mu t- In that model the selective 
advantage a/ (3 of the master sequence is given by the number of offsprings, iV re p, while (3 
is given by the reproduction probability Prep- Since they observed the extinction of the 
master sequence as error threshold, their data must be compared with Eq. @, thus 

-^mut = 1 ~~ 



iV r ep 



corresponds to their error threshold transition and 
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to the mutational meltdown. The phase diagram of this model is reported in Figure 2. The 



numerical data of Ref [I3| are well reproduced, except for a small shift of the error threshold 
transition. 

In summary we find that for a very steep fitness function (a 3> (3), by increasing the 
mutation rate \i or the genome length L one always observes the error threshold; but for 
moderately fitness difference among strains the mutational meltdown causes the extinction 
of the whole population while retaining the quasispecies distribution. 
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FIG. 1. Phase diagram for the error threshold and mutational meltdown transitions. The 
continuous line corresponds to the error threshold q e , Eq. ([?]), the dashed line to the mutational 
meltdown q m , Eq. (S), with = 0.5, the dotted line to the mutational meltdown with f3 = 0.5, the 
dashed-dotted line to the mutational meltdown with (5 = 1.0, corresponding to Muller's ratchet q r , 
Eq. (|9|). Since q = 1 — /xL, the increasing of the mutation rate or the genome length corresponds 
to lowering q. The lower left corner correspond to extinction. 
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FIG. 2. Phase diagram for the error threshold (Muller's ratchet version, vertical lines) and 



mutational meltdown (curves) transitions for model M2 of Ref. 
q = 1 — -P mu t- Continuous line iVrep = 2, dashed line iVrep = 5, dotted line N Te p = 10. 
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